Retardation effects in the Holstein-Hubbard chain at half-filling 
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The ground state phase diagram of the half-filled one-dimensional Holstein-Hubbard model con- 
tains a charge- density- wave (CDW) phase, driven by the electron-phonon (e-ph) coupling, and a 
spin-density-wave (SDW) phase, driven by the on-site electron-electron (e-e) repulsion. Recently, 
the existence of a third phase, which is metallic and lies in a finite region of parameter space 
between these two gapped phases, has been claimed. We study this claim using a renormalization- 
group method for interacting electrons that has been extended to include also e-ph couplings. Our 
method Q treats e-e and e-ph interactions on an equal footing and takes retardation effects fully 
into account. We find a direct transition between the spin- and charge-density wave states. We 
study the effects of retardation, which are particularly important near the transition, and find that 
Umklapp processes at finite frequencies drive the CDW instability close to the transition. We also 
perform determinantal quantum Monte Carlo calculations of correlation functions to confirm our 
results for the phase diagram. 
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The interplay between electron-electron (e-e) and 
electron-phonon (e-ph) interactions leads to important 
effects in low-dimensional materials such as molecular 
crystals, charge transfer solids [2], conducting polymers 
01, and fullerenes Q. In narrow band electronic materi- 
als, perhaps the simplest model capturing this interplay 
is the Holstein-Hubbard model (HHM), where the e-e in- 
teractions are described by a on-site repulsive Coulomb 
term, and the electrons are coupled to dispersionless op- 
tical phonons in localized vibrational modes [3. 

In the one-dimensional HHM (1DHHM) at half-filling, 
early quantum Monte Carlo (QMC) calculations [6] sug- 
gested that there are only two phases: the Peierls charge- 
density-wave (CDW) and the Mott spin-density-wave 
(SDW) state. The boundary between these two phases 
was predicted to lie along the line in parameter space 
where an "effective" e-e interaction vanishes: f7 e ff = 
U — 2g^ p /uJo — 0, where U is the Hubbard on-site e-e re- 
pulsion, g ep is the electron-phonon coupling, and uoq is the 
phonon frequency. More recently, several authors have 
proposed that a third phase might exist near U e E — 0: a 
metallic, Luttinger liquid, phase B H, B , or an off-site 
pairing superconducting phase Large scale QMC 

studies have indicated that there is a metallic region 
with dominant superconducting (SC) pairing correlations 
between the CDW and SDW regions. DMRG studies H3 
suggest that SC does not exist but instead that both the 
spin and charge gaps vanish only for C7 e ff — 0, suggest- 
ing that a metallic phase (with no dominant SC correla- 
tions) may exist only exactly on the boundary between 
the CDW and SDW phases. This is also the conclusion of 
two-step renormalization-group studies [14] and Lanczos 
diagonalization To attempt to determine which of 

these scenarios is correct, we study the problem here us- 
ing a recently developed extended renormalization group 



approach 

At half-filling, Umklapp scattering creates a strong 
tendency to open a charge gap. From the perspective of 
weak-coupling approaches, it is highly non-trivial to have 
a finite metallic, or SC, region. If such a phase is to ex- 
ist, it must be that the dynamical nature of the phonons 
effectively suppresses Umklapp scattering. Therefore, re- 
tardation effects must be taken into account in order to 
investigate this issue. For this purpose, we use a multi- 
scale functional renormalization-group (MFRG) method 
[J. Our MFRG is an extension of the RG for interact- 
ing fermions [15] that are also coupled to bosonic modes 
and applies to both weak (A « 1) and strong (A ^> 1) 
electron-phonon coupling limit (A = 2N(0)g^ p /uJo, N(0) 
is the electron density of states at the Fermi level). For 
a spherical Fermi surface, the MFRG reproduces Eliash- 
berg's theory at the SC instability |lj, and it has also 
been applied in the study of effects of phonons in ladder 
systems |l6j . 

The 1DHHM is given by the Hamiltonian 



H = 



H.c. 



U^Tn 



+#e P + CLijn^a + ^ a\ai, 
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where c\ a (ci l(T ) is an electron creation (annihilation) op- 
erators at site i with spin <j, is the electron number 
operator, a\ (a$) is a creation (annihilation) operator for 
an optical phonon at site z, tis the nearest-neighbor elec- 
tron hopping integral. We use units such that t = 1 = ft. 

Using a path integral formulation and integrating out 
the phonon fields exactly, we find that the effective (re- 
tarded) e-e interaction becomes Q: 
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where k = (fc, uo). We use a notation in which, after scat- 
tering, an incoming electron with momentum and fre- 
quency k x (fc 2 ) goes out with k 4 (fc 3 ), so that k x + fc 2 = 
/c 3 + /c 4 . In the anti-adiabatic limit, where c^o — > oo, all 
the electronic frequency dependences are suppressed, and 
the HHM maps onto the standard Hubbard model with 
a renormalized U e s. At half-filling, its ground state is 
charge-gapped SDW for repulsive interactions and spin- 
gapped degenerate CDW/SC for attractive interactions. 
The transition between SDW and degenerate CDW/SC 
occurs when the bare coupling changes sign, that is when 

U e ft = 0. 

In the MFRG approach at the one-loop level, the RG 
flow equations for the coupling functions, g(k 11 k 2 ,k 3 , k 4 ) 
with initial conditions given by (|2j), are given by [j: 

dgik^k^k^ 



dA 

d 

d P-jX [ g a(p)G a , k 2 ,k)g(p, k, k 3 ) 



-dA 

dp ^ [G A (p) Ga (q 2 )][ -2g (k^p.q^g (q 2 ,k 2 ,k 3 ) 
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where k = k ± + k 2 - p, = P + k 3 - k x . q 2 = p + 
k.3 — k.2i f dp = J dp"}2 UJ l/(27rP) J and Ga is the self- 
energy corrected propagator at energy cut-off A. Since 
the interaction vertices are frequency dependent, there 
are also self-energy corrections. At the one-loop level, 
the self-energy MFRG equation is: 



dA 



dp^ [G A (p)} [2g(p, k, k) - g(k, p,k)}. (4) 



We have solved the coupled integral-differential equa- 
tions, (0) and Q, numerically with two Fermi points 
(iVfc = 2) and by dividing the frequency axis into fifteen 
segments (N u = 15). Fig. [I] shows the discretization 
scheme for Nk = 2 and = 15. 

We next calculate within our MFRG approach the RG 
flow of susceptibilities in the static (zero frequency) and 
long- wavelength limit. In particular, the SC susceptibil- 
ity is given by: x| C (0, 0)=/L>(l, 2)(c Pl , i c_ Pl , T c 



"P2,T P2,| /5 



and the SDW and CDW susceptibilities can be writ- 
ten as: xi(7r,0) =/^(l,2)(ct ijaiCpi ^ jC7l c^ >cr2 Cp 2)0 . 2 >, 
where p$ is the momentum at energy J D(l,2) = 

il^|>A^i J (^)i|^|>A ^2^(6) E^,^ ^i^ 2J and J (0 
is the Jacobian for the coordinate transformation from k 
to f fe . For * = SDW: s T = 1, s x = -1, and for S = CDW: 
= l,Sj = 1. The dominant instability is determined 
by the most divergent susceptibility as the cut-off A is 
lowered. The RG flow for the SC susceptibility is given 
by: 



dA 



--Jdp-^[G A (p)G A (-p)](Zl c (p)f 




Figure 1: Discretization of the momenta in the Brillouin zone 
and frequencies in the frequencies axis. This figure shows the 
case N k = 2, N w = 15. 



dZl c {p)_ r d 
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where g sc (p',p) = g(p',—p',—p), and MFRG flows for 
the SDW and CDW susceptibilities are, 



^Xa(^O) 
dA 



(5) 
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where Q (tt,0). For S SDW: g 6 (p^_,p) = -g(p + 
Q,y,p), and for S = CDW : g 6 (^p) = 2g(^p + Q,p) - 
g(p + Q : p' \p). The function Z 5 (p) is the effective vertex 
in the definition of the susceptibility x 6 ■ Its initial RG 
value is 1. The MFRG equations for susceptibilities are 
solved with initial condition Xa=a = 0- 

In g-ology 0, Eij there are only four couplings, 
corresponding to forward (#2, #4), backward (#i), and 
Umklapp (#3), scattering. The charge and the spin 
parts are governed by gs and gi, respectively. Under 
the MFRG, each one of these couplings carries frequency 
dependence, gi (uu\ , 002 5 ^3) • In the weak e-ph coupling 
limit (A <C 1), the two-step RG is a good approximation, 
and the couplings are separated into two types: high 
frequency transfer, \uo\ — c^4 1 > c^o, and low frequency 
transfer, \uj\ — < coq. However, our MFRG analysis 
reveals that the couplings develop additional non-trivial 
frequency dependence, particularly when the e-ph cou- 
pling is comparable to the e-e coupling and C7 e ff ~ 0. 
As we shall see, understanding this frequency structure 
is critical to resolving the current controversy about the 
behavior in the region near the CDW-SDW transition. 

Deep inside the CDW and SDW regions, we fix ujq = 
1.0 and U = 0.5, and show results of the RG flows for the 
susceptibilities and couplings for different values of g ep . 
For small e-ph coupling (g ep = 0.2, and U e s > 0), the 
SDW susceptibility exhibits a strong divergence, while 
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Figure 2: Left: flows of SC, SDW, and CDW susceptibilities 
for U = 0.5 and ooo = 1.0. Right: flows of Umklapp #3 and 
back-scattering #i, at zero frequencies. Top: g ep = 0.2 (U e s > 
0). Bottom: g ep = 0.8 (J7 e ff < 0). 

both CDW and SC susceptibilities are suppressed (Figj2j 
top). This is expected, since the on-site repulsion dom- 
inates over the retarded attractive interaction mediated 
by the phonons. A charge gap develops, with no spin 
gap, which can be inferred from the flow of the couplings: 
Umklapp (#3) diverges, whereas back-scattering (gi) does 
not. For large e-ph coupling (# ep = 0.8, and f7 e ff < 0), 
the CDW susceptibility diverges (Fig. bottom). Now 
there are both spin and charge gaps, and, correspond- 
ingly, both Umklapp (gs) and back-scattering (gi) are 
divergent. 



plings is qualitatively the same as in the rest of the 
SDW phase (Fig. [2J top). The only difference is that 
the gap decreases and eventually goes to zero at the 
transition. Fig. OH shows the flows for # ep = 0.55 (C/ e ff 
slightly above zero). The SC susceptibility becomes en- 
hanced, but the CDW susceptibility still dominates. In- 
terestingly, #i(0, 0, 0) diverges but #3(0, 0, 0) does not. In 
ID problems without retardation, the usual interpreta- 
tion is that the CDW instability occurs when #1 — > — 00 
and gs — > — 00 0, [HI I^J. In the present case, since 
#3(0,0,0) — > 0, we need to look at the frequency depen- 
dence of the couplings in order to understand what is 
driving the CDW instability. 

In the MFRG approach, we obtain the RG flow of all 
the #2(^1,6^2,^3) couplings and self-energies, and there- 
fore can analyze how this frequency dependence evolves 
with the RG flow. Consider first the cases deep in the 
SDW and CDW phases. Fig. @] shows contour plots of 
#3(^1,^2,^2,^1) which corresponds to an Umklapp pro- 
cess with zero-frequency transfer, \ui — 6^4 1 = 0. We plot 
the value of the coupling at an RG scale i right before 
the critical scale i c when the instability occurs. For the 
SDW phase (Fig. 21 left), the existence of a charge gap is 
signaled by divergence in the Umklapp channel, and the 
most divergent #3 couplings are the ones close to zero 
frequency. Deep inside the CDW phase, #3(0, 0, 0, 0) also 
diverges, as we have seen before from Fig. [21 However, 
the most divergent couplings are for large values of uj\ 
and UO2 (see Fig. @J • 
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Figure 3: Left: flow of susceptibilities for U = 0.5, ooo = 1.0, 
# ep = 0.55 (£4ff < 0). Right: flows of the Umklapp scattering 
#3 and back-scattering #i at zero frequency. 

We next consider the region close to the CDW-SDW 
transition where U e s ~ 0. For U e s slightly below zero 
(# ep = 0.48), the behavior of susceptibilities and cou- 



Figure 4: Plots of the Umklapp scattering #3(^1,^2,^2,^1) 
for U = 0.5, and u = 1.0. Left: # ep = 0.2. Right: # ep = 0.8. 

The situation for # ep = 0.55, shown in Fig. is more 
intriguing. Umklapp scattering is renormalized to large 
values in most part of the frequency space. However, for 
frequencies near zero Umklapp scattering flows to very 
small values. From the RG flow of the susceptibilities 
(Figs. [21 and 0), it is clear that there is CDW instability 
for U e ff > and a direct transition from CDW to SDW. 
From the frequency dependence of #3 we conclude that 
close to the transition to the SDW, the CDW instability 
is being driven by Umklapp processes at high frequencies. 
These are processes at small frequency transfer, \ui — 
ou a I ~ < c^o but that nevertheless involve electrons with 
high frequencies (uj\ and 0J2). In a two-step RG analysis, 
the couplings #3(^1,^2,^2,^1), with different uj\ and 0J2 
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Figure 5: Plot of the Umklapp scattering ^3(^1, UJ2, ^2, 
for U = 0.5, 000 = 1.0, and g 3 = 0.55. Note that # 3 (0,0,0) is 
flowing towards zero. 



are all indistinguishable since \ui — c^4 1 = for all of 
them. Clearly, the two-step RG fails in this region. 

As an independent (partial) confirmation of our MFRG 
results, we have also performed determinantal QMC 
calculations for the Holstein model (U = 0). For the 
charge exponent, Kqbw = lim q ^o 7rS p (q)/q, we obtain 
that i^cDW > 1 when g ep is smaller than some value 
that depends on coo- This result agrees with that ob- 
tained in [llj, using stochastic series expansion QMC 
jijj. For a Luttinger liquid, the scalings of ground state 
correlation functions are determined solely by the charge 
(K p ) and spin (K a ) exponents. For example, in the spin- 
gapped regime, where K a = 0, CDW and SC correla- 
tion functions scale as CDW (x) oc x~ aK ? = x~ KcjDW , 



and O nj) oc x 



-P/k p 



" Ksc , with a 



1 



IlllSEi. The dominant correlation is of CDW (SC) 
type for K p < 1 {K p > 1). This relation is not guaran- 
teed to hold in the presence of phonons and retardation 
effects (Hj. 
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Figure 6: SC and CDW correlations for 38-sites Holstein 
model (cj = 1.0, g ep = 0.5), with i^cDw = 1.032 ± 0.005. 

Using the determinantal QMC allows us to calculate 
the pairing and charge correlations directly (Fig. EJ). 
We find that the charge correlation function decays more 



slowly. This provides, at least for the case U = 0, confir- 
mation of our MFRG results and strongly suggests that 
there is no region of dominant SC correlations in the half- 
filled 1DHHM, even though the scaling exponent of the 
charge correlation function can be larger than 1. 

In conclusion, we have studied the ground state of 
1DHHM at half-filling using the MFRG method. This 
technique enables us to treat retardation effects from the 
phonons in a systematic way. We find SDW and CDW 
phases, and a direct transition between them. Analysis 
of the frequency dependence of the gs shows a shift in 
spectral weight indicating that the CDW instability near 
the transition is driven by dynamical Umklapp processes. 
Our determinantal QMC results for the charge exponent 
and correlation functions for the Holstein model confirm 
our MFRG predictions and suggest that having a charge 
exponent larger than one for finite size system does not 
mean dominant SC correlations because of breakdown of 
TLL relations due to retardation. 
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